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(57) Abstract 

The invention relates in general to digital encoding and decoding of information. More particularly, the invention relates 
to efficient implementation of digital analysis and synthesis filter banks used in digital encoding and decoding. In a preferred 
embodiment of the invention, the length of the filter bank used to implement critically-sampled analysis and synthesis filter 
banks may be adaptively selected. 
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DESaOFnON 

LOW COMPUTATTONAL-COMPLEXITY DIGITAL FILTER BANK 

5 

> 

Technical Field 

The inveiition relates in general to digital encoding and decoding of information. More 
10 particularly, the invention relates to efficient implementation of digital analysis and synthesis filter 
banks used in digital encoding and decoding. In a preferred embodiment of the invention, the length 
of the filter bank used to implement critically-sampled analysis and synthesis filter banks may be 
ad^tively selected. 

Throughout the following discussion and especiaUy in the background discussion, more particular 
15 mention will be made of audio applications; however, it should be understood that the present 

invention is applicable to a range of applications wider than just that of audio encoding and decoding. 

Background Art 
Introduction 

20 There is considerable interest among those in the field of signal processing to develop efficient 
means to transmit or store information. Improving coding efficiency includes (1) reducing 
informational requirements, that is, reducing the amount of information required to adequately 
represent a signal during transmission or storage, and (2) reducing processing requirements, that is, 
reducing the amount of processing required to implement the encoding and decodmg processes. 

25 In high-quality audio coding applications, informational requirements can sometimes be reduced 
without loss of perceptible audio quality by exploiting various psychoacoustic effects. Signal 
recording, transmitting, or r^roducing techniques which divide the useful signal bandwidth into 
narrow bands with bandwidths approximating the human ear's critical bands can exploit 
psychoacoustic masking effects. Such techniques divide the signal bandwidth with an analysis filter 

30 bank, process the signal passed by each filter band, and reconstruct a replica of the original signal 
with a synthesis filter bank. 

Two common coding techniques are subband coding and transform coding. Subband coders and 
transform coders can reduce the informational requirements in particular frequency bands where the 
p noise caused by the resulting coding inaccuracy is psychoacoustically masked. Subband coders may 

35 be implemented by a bank of digital bandpass filters defining subbands of varying bandwidth. 
Transform coders may be implemented by any of several time-domain to frequency-domain 
transforms. One or more adjacent transform coefficients are grouped together to define "subbands" 
having effective bandwidths which are sums of individual transform coefficient bandwidths. 

The mathematical basis for digital subband filter banks and digital block transforms is essentiaUy 
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the same. See Tribolet and Ciochiere, "Frequency Domain Coding of Speecli," l^^tf Trans. 
AcQust.. Speech, and Signal Proc., ASSP-27, October, 1979, pp. 512-30. Therefore, throughout the 
following discussion the concq)ts associated wifli terms such as "subband coder" and "transform 
coder" gOTerally apply to both a true subband coder and a transform coder. The term "subband" 
5 refers to portions of the useful signal bandwidth ^^er implemented by a true subband coder or a 
transform coder. The terms "transform" and "transforming" include digital filters and digital 
filtering, respectively. 

In most digital coding ai5)lications, processing requirements can be reduced by increasing the 
efficiency of subband filtering. Improved processing efficiency permits implementation of encoders 
10 and decoders \^*ich are less expensive to build, or vMch impose lower signal propagation delays 
through an encoder/decoder system. 

In many transform coder systems, the analysis and synthesis filter banks are implemented by 
discrete time-domain to frequency-domain transforms such as the Discrete Fourier Transform (DFT), 
the Discrete Cosine Transform (DCT), and the Discrete Sine Transform (DST). The number of time- 
15 domain signal samples, referred to herein as the time-domain signal sample block length, processed 
by such transforms is sometimes called the transform lengtii, and the amount of processing required to 
perform these transforms is generally proportional to the square of the time-domain signal sample 
block length. 

The number of frequency-domain transform coefficients generated by a transform is also 

20 sometimes called the transform length. It is common for the number of frequency-domain transform 
coefficients gaierated by the transform to be equal to the time-domain signal sample block length, but 
this equality is not necessary. For example, one transform referred to herein as the E-TDAC 
transform is sometimes described in the art as a transform of length N that transforms signal sample 
blocks with a lengfli of 2iV samples. It is possible, however, to also describe the transform as one of 

25 length \^4iich generates only unique ftequency-domain transform coefficients. Thus, in this 
discussion the time-domain signal sample block length and the discrete transform length are generally 
assumed to be synonyms. 

Various techniques have been utilized to reduce the amount of time required to perform a 
transform, or to reduce the processing power required to perform a transform in given amount of 

30 time, or both. One technique is taught in Narasimha and Peterson, "On the Computation of the 
Discrete Cosine Transform," IEEE Trans, on Communications , COM-26, June, 1978, pp. 934-36. 
Briefly, this technique evaluates an ^-point DCT by rearranging or "shuffling" the samples 
r^resenting the input signal, performing an ^-point DFT on the shuffled samples, and multiplying 
the result with a complex function. It is ^proximately twice as efficient as other techniques using a 

35 2^-point FFT; however, Narasimha and Peterson only teach how to improve Uie efficiency of filter 
banks implemented by one particular DCT. 

Another technique \^ch yields proximately a two-fold increase in processing efficiency 
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concurrently performs two real- valued discrete transforms of length N with a single complex-valued 
FFT of length N. A transform coder utilizing this technique to concurrently perform a modified DCT 
with a modified DST is described in International Patent Application PCT/US 91/02512, Publication 
No. WO 91/16769 (published October 31, 1991). The significance of these particular modified DCT 
S and modified DST is discussed in Princen and Bradley, "Analysis/Synthesis Filter Bank Design Based 
on Tmie Domain Aliasing Cancellation," EEEE Trans, on Acoust- Speech, Signal Proc . ASSP-34, 
1986, pp. 1153-1161. The authors describe a specific ^plication of these transforms as the time- 
domain equivalent of an evenly-stacked critically-sampled single-sideband analysis-synthesis system. 
They are referred to collectively herein as the Evenly-stacked Time-Domain Aliasing Cancellation 
10 (E-TDAC) transform. 

Another technique to reduce processing requirements is tau^t by Malvar, "L^ped Transforms for 
Efficient Transfonn/Subband Coding," IEEE Trans. Acoust., Speech. Signal Proc . ASSP-38, June, 
1980, pp. 969-78. This technique implements an ^-point modified DCT by performing a V&iV-point 
DST after combining pairs of the samples representing the input signal, or "folding" the N input 
15 signal samples mto a smaller set of %iV points* It is ^proximately twice as efficient as performing 
the modified DCT in a straight-forward manner; however, Malvar only teaches how to fold input 
samples for a filter bank implemented by one specific modified DCT whose input samples have been 
weighted by a specific sine-tapered analysis window. 

The specific modified DCT implemented by Malvar is discussed in greater detail by Princen, 
20 Johnson, and Bradley, "Subband/Transform Coding Using Filter Bank Designs Based on Time 
Domain Aliasing CanceUation," ICASSP 1987 Conf. Proc. , May 1987, pp. 2161-64. The authors 
describe this transform as the time-domain equivalent of an oddly-stacked critically sampled single- 
sideband analysis-synthesis system. It is referred to herein as the Oddly-stacked Time-Domain 
Aliasing Cancellation (O-TDAC) transform. 
25 It is desirable to implement encoders and decoders with the ability to use different time-domain 
signal sample block lengths in order to optimize coder performance. It is well known in the art that 
longer time-domain signal sample block lengths improve the selectivity or frequency-resolving power 
of transform coders, and better filter selectivity generally improves the ability of a transform coder to 
exploit psychoacoustic masking effects. 
30 But longer time-domain signal sample block lengths degrade the time-resolution of a subband filter 
bank. Inadequate time-resolution can produce audible distortion artifacts when quantizing errors of 
signal events, such as transients, producing pre-transient and post-transient ringing which exceed the 
ear's temporal psychoacoustic masking interval. Hence, it is desirable that techniques which improve 
subband filter bank processing efficiency should also permit adaptive selection of the time-domain 
35 signal sample block length. 

The importance of time-domain signal sample block length and its effect upon filter bank 
frequency-domain resolution and time-domain resolution is discussed in more detail in International 



wo 92/22137 I'd/ us^vz/iw/o / 

-4- 

PMent Application PCT/US 91/02512, Publication No, WO 91/16769, cited above, which is hereby 
incozporated by reference in its ratirety. 



5 ^^ffrJfffrV? nf Invention 

It is an object of the present invention to provide for a subband/transfonn encoder and a 
subband/tiansform decoder of digital information by means of analysis filtering and synthesis filtering 
requiting lower processing requirements, or imposing lower processing delays, or both. 

It is another object of the present invention to provide for a subband/transform encoder and a 

10 subband/transfonn decoder of digital information requiring lower processing requirements, or 
imposing lower processing delays, or both, by means of analysis filtering and synthesis filtering 
permitting ad^tive selection of the filter-bank length. 

Further details of the above objects and still other objects of the invention are set forth throughout 
this document, particulariy in the Modes for Carrying Out the Invention, below. Although the 

IS invention is more particulady described for audio encoding and decoding plications, it should be 
appreciated that the invention is much broader and may be implied to other applications. Throughout 
this Description, discussion of encoders incorporating the present invention also pertains more 
generaUy to signal-analysis filtering plications, and discussion of decoders incorporating the present 
invention pertains more graerally to signal-synthesis filtering ^plications. 

20 In accordance with the teachings of the present invention in one embodiment, an encoder provides 
for the encoding of input signal samples representing a time-domain signal. The input samples vMch 
are weighted by an analysis^window function are buffered into time-domain signal sample blocks. 
Pairs of signal samples in the time-domain signal sample blocks are combined by a forward pre- 
transform function to generate modified samples. Frequency-domain transform coefficients are 

25 generated by applying a discrete digital transform to the modified samples. Spectral information is 
generated by applying a forward post-transform function to the frequency-domain transform 
coefficirats. 

Also in accordance with the teachings of the present invention in one embodiment, a decoder 
provides for the decoding of digitally encoded spectral information. Frequency-domain transform 

30 coefficients are generated by applying an inverse pre-transform function to the spectral information. 
Time-domain transform coefficients are generated by plying an inverse discrete digital transform to 
the firequency-domain transform coefficients. Time-domain signal sample blocks are generated by 
applying an inverse post-transform function to the time-domain transform coefficients, and ou^ut 
samples v^ch correspond to the input samples to a companion encoder are generated by overlping 

35 and adding samples in adjacent time-domain signal sample blocks. 

The various features of the present invention and its preferred embodiments are set forth in greater 
detail in the following Modes for Carrying Out the Invention and in the accompanying drawings. 
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Brief Description of Drawings 
Figure 1 is a functional block diagram Ulustrating the basic functional structure of an encoder 
incorporating a preferred embodiment of the present invention. 
S Figure 2 is a functional block diagram illustrating the basic functional structure of a decoder 
incorporating a preferred embodiment of the present invention. 

Figure 3 is a flowgr^h illustrating the forward pre-transform function applied to a 16-sample 
time-domain signal sample block to form an 8-sample modified sample block for a basic embodiment 
of the present invention permitting implementation of an E-TDAC transform analysis filter bank by a 
10 DCTandDST, 

Hgure 4 is a flowgr^h illustrating the forward pre-transform function applied to a 16-sample 
time-domain signal sample block to form an 8-sample modified sample block for an alternative 
embodiment of the present invention permitting implementation of an E-TDAC transform analysis 
filter bank by a DFT. 

IS Figure S is a hypothetical gr^hical r^resentation illustrating the time-reversal regions of the 
time-domain aliasing component created by the E-TDAC transform using the conventional TDAC 
phase term. 

Hgure 6 is a hypothetical graphical representation illustrating the time-reversal regions of the 
time-domain aliasing component created by tiie E-TDAC transform using the TDAC phase term 
20 required to cancel time-domain aliasing in an A^sample length block overl^ped with a subsequent 
V^ys^sample length block. 

Figure 7 is a hypothetical gn^hical representation illustrating the boundary between time-reversal 
regions of the time-domain aliasing component in a ViiV-sample length block. 

Figure 8 is a hypothetical gnq}hicai representation of a bridge transform illustrating the time- 
25 reversal regions of the time-domain aliasing component. 

Figure 9 is a flowgraph illustrating the forward pre-transform function applied to a 16-sample 
time-domain signal sample block to form an 8-sample modified sample block permitting 
implementation of an adaptive-lengUi E-TDAC transform analysis filter bank by a DCT and DST. 
Figure 10 is a flowgrqih illustrating the forward pre-transform function sailed to a 16-sample 
30 time-domain signal sample block to form an 8-sample modified sample block permitting 
implementation of an adaptive-lengtii 0-TDAC transform analysis filter bank by a DST. 

Modes for Carrying Out the Invention 
35 L Overview of Functional Structure 

Figure 1 illustrates the basic functional structure of a transform-based encoder incorporating an 
embodiment of the present invention. According to this structure, an encoder comprises buffer 102 



wo 92/22137 rc i / uay^/ u*# / o / 

-6- 

vMch buffers input samples received from input path 100 into time^iomain signal sample blocks, 
forward pre-transform 106 ^ch generates modified samples by combining pairs of signal samples 
received ftom buffer 102 and in response to information received firom path 104 establishing the 
number of signal samples constitutmg a time-domain signal sample block, forward transform 108 
5 vMch transforms the modified samples into frequency-domain transform coefficients by applying a 
transform vfhose length is ad^ted in response to information received from path 104, forward post- 
transform 110 vMch generates spectral information from the frequency-domain transform coefficients 
and in response to the information received from path 104, and formatter 112 which assembles digital 
information including the spectral information into a form suitable for transmission or storage along 
10 path 114. The functions performed by buffer 102 and formatter 112 are not discussed in detail 
herein. 

Figure 2 illustrates the basic functional structure of a transform-based decoder incorporatmg an 
CTibodiment of the present invention. According to this structure, a decoder comprises deformatter 
202 vMcb extracts spectral information and information establishing the inverse transform length 

15 from the encoded digital signal received from path 200, inverse pre-transform 206 which genwates 
frequency-domain transform coefficients from the extracted spectral information and in response to 
the information establishing the inverse transform length received along path 204, inverse transform 
208 v^iiich transforms the frequency-domain transform coefficients into time-domain transform 
coefficients by flying a transform whose length is ad^ted in response to information received from 

20 path 204, inverse post-transform 210 which generates signal samples from the time-domain transform 
coefficients and in response to information received from path 204, and output processor 212 which 
generates along path 214 ou^t samples corresponding to the input samples to a companion encoder 
in response to the signal samples. The fonctions performed by deformatter 202 and output processor 
212 are not discussed in detail herein. 

25 It should be ^rpreciated from a study of the foUowing disclosure and the accompanying claims that 
some elements shown in Figures 1 and 2 are not required to practice various embodiments of the 
present invention. 

A basic embodiment of the present invention is introduced in some detail before alternative 
embodiments are discussed. This basic embodiment uses fixed-length E-TDAC transforms to 
30 implement the analysis and synthesis filter banks. Preferred embodiments of various features are 
described throughout the discussion. 

n. Basic Embodiment of Invention 
A. Input Sample Buffering 
35 A buffer, rqjresented by box 102 in Figure 1, receives signal samples and groups them into a 
sequence of time-domain signal sample blocks. Each block comprises signal samples. The signal 
samples may be received from the sampling of an analog signal, from the generation of samples 
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lepresenting or simulating an analog signal ^ or from any other source of discrete- valued samples 
which correspond to a time-domain signal. 

It is well known in the art that the frequency-resolving power or selectivity of a filter bank 
implemented by a discrete transform improves as the transform length increases. It is also well 
5 known that filter selectivity may be affected significantly by weighting the time-domain signal samples 
by a weighting function commonly called a window. See generally, Harris, "On the Use of Windows 
for Harmonic Analysis with the Discrete Fourier Transform," Proc. IEEE , vol. 66, January, 1978, 
pp. 51-83. 

Hie E-TDAC transform used in the basic embodiment of the present invention requires window 
10 weighting, both weighting of the time-domain signal samples in an encoder prior to forward transform 
filtering, referred to as analysis windowing, and weighting of the recovered time-domain signal 
samples in a decoder after inverse transform filtering, referred to as synthesis windowing. Analysis- 
and synthesis-window weighting are discussed below only briefly. It is assumed herein that the 
buffered signal samples are weighted by an analysis window as may be required or desired. Input 
15 signal samples may be weighted by an analysis window prior to or subsequent to their receipt by the 
buffer without dq^arting from the scope of the present invention. 



B. Analysis Filter Bank - Forward Transform 
Although the forward pre-transform function discussed below is applied to time-domain signal 
20 sample blocks prior to application of the forward transform, it is necessary to introduce the forward 
transform before the forward pre-transform function can be fully described. The forward transform is 
rqjresented by box 108 in Figure 1, 

The E-TDAC transform used in the basic embodiment of the present invention is equivalent to the 
alternate application of a Modified Discrete Cosine Transform (MDCT) with a Modified Discrete Sine 
25 Transform (MDST). The MDCT and the MDST, shown in equations 1 and 2 respectively, are 

C(k) = J^jc(ii)cos(2TJfcii^) forO ^ A < TV (1) 

Sik) - Tx(n)sin(2Tk!Ll^) forO ^ * < TV (2) 

ii-o TV 

where k = frequency-domain transform coefficient number, 
n = time-domain signal sample number, 
30 TV = time-domain signal sample block length, 

m = phase term required for TDAC (see equation 6), 
x(n) — time-domain signal sample n, 

C(k) = MDCT frequency-domain transform coefficient ky and 
S(k) = MDST frequency-domain transform coefGcient k. 
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The E-TDAC transform produces one of two alternating sets of ftequency-domain transform 
coefBcients in response to each time-domain signal sample block. These sets of frequency-domain 
transform coefiGcients are of the form 



{cm, 



Qk) for 0 ^ it < ^ 

0 imk^ — 
2 



(3) 



S(k) forO < k ^ ^ (4) 
0 for ifc » 0 

v/h&co i = time-domain signal sample block number. Each set of coefBcients generated by the 
MDCT and the MDST are referred to herem as MDCT coefficient sets and MDST coefficient sets, 
respectively. 

Princen and Bradley showed that with the proper phase term m and a suitable pair of analysis- 
10 synthesis windows, the E-TDAC technique can accurately recover an input signal from an alternating 
sequence of overi^qped fixed-length MDCT coefficient sets and MDST coefficient sets of the form 
{C(fc)}o, {5(ft)}t, {C(*)h, {S(k)}^. . • . - 

Using only the alternate MDCT coefficient sets and MDST coefficient sets produces a time- 
domain aliasing component^ but the aliasing component may be cancelled by choosing the appropriate 

15 lAase term m for equations 1 and 2, applying the forward transform to overtyped analysis-window 
weighted time-domain signal sample blocks, and by synthesis-window weighting and adding adjacent 
overlq>ped time-domain signal sample blocks recovered by the inverse transform. 

The phase term m in equations 1 and 2 controls the phase shift of the time-domain aliasing 
distortion. To cancel this alias distortion and accurately recover the original time-domain signal, 

20 E-TDAC requires the aliasing to be as follows: for the MDCT, the time-domain alias component 
consists of the first half of the sampled and windowed signal reversed in time about the one-quarter 
point of the sample block and the second half of the sampled and windowed signal reversed in time 
about the three-quarter point of the sample block; for the MDST, the alias component is similar to 
that for the MDCT except its amplitude is inverted in sign. These relationships arc illustrated in 

25 Figure 5 in \siiich the time-domain aliasing component, shown by a broken line, and the desired 
signal, shown by a solid line, have been weighted by a synthesis window. 

The i^iase term required to produce the appropriate aliasing components for alias cancellation is 
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C. Forward Pre-Tiansform Function 
The pnx:essing requirements of the technique used to evaluate the MDCT and the MDST may be 
reduced by ^plying a forward pre-transfonn function to the time-domain signal sample bloclcs to 
produce blocks of modified samples, and applying a V^Mpoint DCT and a V^^-point DST to the 
5 modified sample blocks for the iV-point MDCT and MDST, respectively. The forward pre-transform 
fimction is r^resented by box 106 in Figure 1. For E-TDAC, the pre-transform fimction combines 
pairs of signal samples in each fime-domain signal sample block of length N to produce a block of 
modified samples of length V^^. 

The mathematical basis for using a ViiV-point DCT applied to modified samples to perform the 
10 iV-point MDCT may be seen by first substituting equation 6 for the phase term m into equation 1. 
The MDCT in equation 1 may be expr^sed as 

5:*(«)cos(i^I«+^+^]) + (7b) 
"'t 



£ jc(fl)cos(i^In+^+^]) + (7c) 



15 £ Jt(n)cos(^ [11+^+4 1) fbrOSk<N. (7d) 

„3w N 4 2 

By setting d = and substituting it into expression 7a, by setting d = — -1-n and substituting 

4 37V ^ 

it into expressions 7b and 7c, and by setting d = n-—. and substituting it into expression 7d, it may 

be seen that equation 1 can be rewritten as 



dk) = £ {x(^+rf) +;c(^-l-^} cos(2^y4]) + (8a) 
d-0 4 4 Jy 2 



20 5^ W^'-t) cos(^M+1]) forO ^ k < M (8b) 

. 4 4 rt i 



"^7 



FinaUy, by defining a new sequence 
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y(n) = l|x(I^+«] mod AO + x([^-l-n] mod AO j for 0 i « < 'AN, 

where [Q mod Af tepteseats the value of i modulo M, the expressions 8a and 8b may be combined 
and written as 



(9) 



C(k) = J^y(/i)cos(^[n+l]) forO ^ < iV (10) 
■•0 N 2 

5 wMchisa%iV-pomtDCTforyf>l'. 

Fiom a similar derivation, it can be shown that the MOST of length can be implemented by a 

DST of length ^AN; 

S(k) - '£z(n)^i^[n^h) forO !S Ar < (H) 

Where z(n) = mod iV) - x([^'l-n] mod Aoj forO ^ « < %M (12) 

10 It should be appreciated that the foiwaid pie-transform function for this basic embodiment, as well 
as for alternative embodiments discussed below, can be performed by any of several implementations 
including software-controlled processors and circuits cqiable of combining pairs of time-domain 
signal samples to form modified samples. One flowgraph for a 16-sample block is shown in Figure 3 
v^ch illustrates the forward pre-transform functions of equations 9 and 12 for a basic embodiment of 

IS the present invention. The minus signs shown within parenthesis denote terms which are 

subtractively combined with an associated sample for the function shown above in equation 12. This 
subtractive combination may be accomplished in circuitry, for example, by negating the value of the 
signal sample rqyresoitations corresponding to the nodes in Figure 3 with minus signs in parenthesis 
and additively combining the resultant representations. 

20 

D. Forward Post-Transform Function 
In principle, the forward transform in the basic embodiment of the present generates frequency- 
domain transform coefficients in response to an input signal viiich are equivalent to the coefficients 
generated by an E-TDAC transform applied to the same input signal. Some alternative embodiments 
25 of the present invention described below require ^plication of a forward post-transform function to 
the coefficirats generated by the forward transform in order to obtain spectral information equivalent 
to transform coefficients generated by a corresponding TDAC transform. 

If an application does not require spectral information, then an encoder incorporating any 
embodiment of the present invention need not s^ply a forward post-transform function to the 
30 frequency-domain transform coefficients generated by the forward transform. For example, the 
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frequency-domaiB tiansform coefficients themselves may be directly transmitted or stored and 
subsequently transmitted to a corresponding receiver for decoding. 

For many ^plications, however, spectral infonnation is required. For example, encoder/decoder 
systems which exploit psychoacoustic principles to reduce coded signal information requirements 
5 usuaUy require spectral information in order to estimate psychoacoustic masking effects of a signal's 
spectral components. 

Frequency-domain transform coefficients generated by the forward transform and spectral 
infonnation generated by the various forward post-transform functions described below are generally 
not suitable for low bit-rate transmission or efficient storage. Various quantization techniques may be 
10 used to reduce informational requirements by taking advantage of a signal's irrelevancy. 

In a practical implementation of an encoder incorporating a basic embodiment of the present 
invention, the forward post-transform function represented by box 110 in Figure 1 may comprise 
quantizing the frequency-domain transform coefficients generated by the forward transform; however, 
quantizing is not required to practice the present invention. 

15 

E. Ou^ut Formatting 

Ou^ut formatting represented by box 112 in Figure 1 is not required to practice the present 
invention but is often used in signal encoding applications. Generally, output formatting assembles 
the spectral information and other information required for transmission or storage. Any additional 
20 side-information needed by a decoder is also assembled into the formatted signal. Frame 

synchronization bits and error detection/correction codes may be used as needed for transmission. 
Database pointers or keys may be added as needed for storage. The formatted data is ready for 
transmission or for storage along path 1 14 shown in Figure 1 . 

25 F. Input Deformatting 

Input deformatting represented by box 202 in Figure 2 is not required to practice the present 
invention but is often used in signal decoding applications. Deformatting extracts spectral information 
and any side information from a formatted signal received from path 200 either by rweipt of a 
transmitted signal or retrieved from storage. 

30 

G. Inverse Pre-Transform Function 
The inverse pre-transform function, represented by box 206 in Figure 2, obtains frequency-domain 
transform coefficients from the spectral information in the received signal. If the spectral information 
in the received signal substantially corresponds to the frequency-domain transform coefficients 
35 generated by an E-TDAC transform, then the inverse pre-transform function in a basic embodiment of 
the present function may be a trivial or essentially null function such as, for example, grouping 
spectral infonnation into blocks. 



wo 92/22137 I'Ul/usyz/iw/o/ 

-12- 

In a practical implementation of a decoder incorporating a basic embodiment of the present 
invention, the inverse pre-transform function may comprise d e qua nt izi n g the encoded digital 
information into a form suitable for input to the inverse transform filter bank; however, dequantizing 
is not required to practice the present invention. 

5 

H. Synthesis Filter Bank - Inverse Transform 
Box 208 in Rgure 2 represents a bank of synthesis filters wliich transforms each set of firequency- 
domain transform coefficients into time-domain transform coefficients. A transform inverse to that 
used in analysis filter bank 108 in Kgure 1 implements synthesis filter bank 208. The inverse 
10 discrete transforms for E-TDAC used in the basic embodiment of the present invention is an 

alternating application of an Inverse Modified Discrete Cosine Transform (IMDCT) and an Inverse 
Modified Discr^ Sine Transform (IMDST) shown m equations 13 and 14, respectively; 

m - - £ C(ife)cos(2xifei!l^) forO ^n<N (13) 



i(rt) « 1 yj^(it)sm(2xfc^) fotO ^n<N 
N to N 



(14) 



15 \^ere C(k) = recovered MDCT frequency-domain transform coefficient 
S(k) = recovered MDST firequency-domain transform coefficient fc, and 
x(n) = recovered time-domain signal sample n. 

L Inverse Post-Transform Function 
20 The processmg requhements of the technique used to evaluate tiie IMDCT and the IMDST may be 
reduced by instead evaluating an Inverse DCT (IDCT) and an Inverse DST (IDST) and applying an 
inverse post-transform function after plication of tiie mverse transforms. This inverse post- 
transform function is rqn:esented by box 210 in Figure 2. 

For E-TDAC, the inverse post-transform function splits time-domain transform coefficients into 
25 signal samples. Using a derivation similar to that discussed above for the forward transform, it can 
be shown that, with an appropriate inverse post-transform function discussed below, the IMDCT of 
length iV can be implemented by an IDCT of lengtii ViiV; 

S(n) = 1 T;a(Jt)C(it)cos(^[/i+l]) forO ^ n < ^AN (15) 
N to N 2 

\siiere ^(n) = recovered time-domain transform coefficient n, and 

N 

2 otherwise. 



1 forit«0,.^ 
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Recovered time-domain signal samples x(n) may be obtained from the time-domain transfonn 
coefBcients Sff^) according to 

je(n) = Hin-^] mod N) for 0 ^ /i < WAT, ^ n < N, (16) 
4 

je(n) = 5>([±I-l-/i] mod AO for UN ^ n < %M (17) 

5 IMth an appropriate inverse post-transform function, the IMDST of length AT can be implemented 
by an IDST of length %7V; 

N 

^^"> " 4 i^ot(k)S(k)sm(2:^[n->-h) forO ^ n < V4A^ (18) 
yv iv 2 

where z(n) = recovered time-domain transform coefficient n. 

Recovered time-domain signal samples may be obtained from the time-domain transform 
10 coefficients z(n) according to 

£(n) = f([/i-^] mod AO forO ^ /i < UN, ^N ^ n < N, (19) 
= -f ([~-l-/t] mod AO for UN ^ n < %N. (20) 



J. Output Sample Processing 

15 An overlap-add process is required by the TDAC transforms to generate samples con^ponding to 
signal samples encoded by a companion encoder. This process, represented by box 212 in Figure 2, 
overl^s adjacent blocks of recovered time-domain samples and adds the samples in one block to 
samples in the adjacent overiapped block. 

The E-TDAC transform used in the basic embodiment of the present invention also requires prior 

20 to overl^add the application of a synthesis window to the recovered time-domain sample blocks. 
The constraints the E-TDAC transform places upon the design of the synthesis window, the analysis 
window, and the ovcriap-add process is discussed fully in the paper by Princen and Bradley referred 
to above. 

25 

m. Alternative ¥lxed-Length Embodiments 
Alternative embodiments of the present invention may achieve greater reductions in processing 
requirements. The following description discusses the differences between these alternative 
embodiments and the basic embodiment described above. 
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A. E-TDAC Implemented by DFT 
In one alternative embodiment of the present invention for an encoder, the forward E-TDAC 
transform is implemented by a Discrete Fourier Transform (DFT). 

A forward pre-transform function generates an alternating sequence of two types of blocks 
5 comprising modified samples; one block type comprising modified samples p(n) and a second block 
type comprising modified samples r(n). Each modified sample is formed from tiie combination of 
OTO pair of signal samples x(n) according to 

p{n) = 1 jx([H+2n] mod N) ^xd^-l -2n] modiV)| for 0 ^ n < ^AN, (21) 

Tin) = ^ jx([^+2n] mod AO -xa^-l-2n] modiV)| for 0 ^ n < (22) 

10 A flowgrsqph for a 16-sample block illustrating fliis forward pre-transform fiinction is shown in 
Figure 4. 

Ite forward E-TDAC transform is implemented by a DFT which generates alternating sets of 
complex-valued fi^quency-domain transform coefBcients P(k) of tiie form T(k)'k'yU(k) and R(k) of tiie 
form V(k)^yW(k) in response to tiie alternating sequence of modified sample blocks; 

15 m = i^p(/i)«''^'*'^ forO ^ Jfc ^ ^AN, (23) 

= E ^^)^^^^''^ for 0 ^ * ^ Vyv (24) 
where j = V-1. 

Spectral information corresponding to E-TDAC transform coefficients C(k) and S(k) is obtained by 
applying a forward post-transform fimction according to; 

20 C(fc) - cosA-m + sin(^) •l^(*), (25) 
N N 

S{k) « sin(lj) • Vik) - cos A • W(Jc). (26) 
N N 

In one alternative embodiment of the present invention for a decoder, the inverse E-TDAC 
transform is implemented by an Inverse DFT (IDFT). 

An inverse pre-transform fimction recovers spectral information C(k) and S(k) corresponding to 
25 E-TDAC transform coefBcients C(k) and S(k), respectively, from tiie encoded signal and generates in 
response to the recovered spectral information an alternating sequence of two types of blocks 
comprising recovered frequency-domain transform coefBcients; one block type comprises recovered 
complex-valued coefBcients P(k) of tiie form 7(i:)+j-t/(*) and a second block type comprises 
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recovered complex-valued coefficients R{k) of the fonn The real and imaginary parts 

of the frequency-domain transform coefficients are obtained according to 

tik) = cos(^) • d(k) + sin(^) • d(!^-k), (27) 
m = sin(^) • C(k) ' cos(^) • C(:^-Jfc), (28) 
5 V{k) = cos(:^)-5(^-*) ^ sin(:^)-5(Jt), (29) 

]«^(*) = sin(-^)-5(:^-*) - cos(:^).^(ifc). (30) 

An inverse transform generates an alternating sequence of two types of blocks comprising 
recovered time-domain transform coefficients by applying an IDFT to the alternating sequence of 
frequency-domain transform coefficient blocks; one block type comprises recovered time-domain 
10 transform coefficients p(k) and a second block type comprises recovered time-domain transform 
coefficients r(k). The IDFT used to recover the time-domain transform coefficients is shown in 
equations 31 and 32; 

Pin) = ^ r Hk)e ^ forO ^ /!< ViiV, (31) 



^'') = i forO ^ /!< %M (32) 

15 Recovered time-domain signal samples i(n) are obtained by applying an inverse post-transform 
function to the alternating sequence of blocks comprising recovered time-domain transform 
coefficients. Signal samples are obtained from blocks comprising the pflc) coefficients according to 

i(n) = Pa—^] mod[^]) = ^([i!!J^] mod[:^]) for 0 ^ « < iV, « even, (33) 

m = ^([-1- ] mod[:^]) = PilEL:^] mod[^]) tot0^n<N,n odd. (34) 

20 Signal samples are obtained from blocks comprising the r(k) coefficients according to 

f(n) = HI—JL] mod[!^]} = mod[:^]) for 0 S « < AT, n even, (35) 
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. -fa-i- ] mod[^]) = -f([3A^^-^»3 mod[^]) for 0 ^ « < A^. u odd. (36) 

2 2 o 2 

B, E-TDAC Implemented by Concurrent DFT 
In another embodiment of the present invention for an encoder, the MDCT and the MDST of one 
5 or more forward E-TDAC transforms are implemented concurrently by one or more DFTs. In single 
channel encoder applications, two adjacent frequency-domain coefficient sets as illustrated in 
expression S above may be generated concurrentiy by a single DFT. In two channel applications, a 
MDCT coefBcient set for channel one may be generated concurrentiy with a MDST coefficient set for 
channel two, immediately followed by a MDST coefBcient set for channel one generated concurrentiy 
10 with a MDCT coefBcient set for channel two. Otiier combinations of coefBcient sets for concurrent 
processing are possible. For additional detail on concurrent transforms, see generally, Brigham, The 
Fast Fourier Transform ^ Englewood Cliffs, NJ: Prentice-Hall, Inc., 1974, pp. 166-67. 

A forward pre-transform function generates a sequence of blocks comprising complex-valued 
modified samples q(n) of tiie form pfnj-^yrfn) M*ere p(n) and r(n) are formed from tiie plication of 
15 the forward pre-transfbrm function described above and shown in equations 21 and 22. 

The MDCT and the MDST constituting tiie forward E-TDAC transform are concurrentiy 
implrairated by a DFT vMch generates compl«-valued frequency-domain transform coefficients Q(k) 
of the form G(k)+yH(k) according to 

= Y^9(n)e ^ forO ^k^ ^AN, (37) 

ii«0 

20 Spectral mformation corresponding to E-TDAC transform coefficients C(k) and S^c) is obtained by 
applying the forward post-transform functions according to 

Qk) . ^|cos(^)-[G(fe) + Gij'k)] * sin(:^)-[H(fc) - |, (38) 

Sik) = ^ jcos(^)-[G(fc) - Gil^'k)] + sm(J^)^[H(k) + /f(^-*)] |. (39) 

In another embodiment of the present invention for a decoder, the IMDCT and the IMDST of one 
25 or more inverse E-TDAC transforms are implemented concurrentiy by one or more IDFTs. 

An inverse pre-transform function recovers spectral information C(k) and S(k) corresponding to 
E-TDAC transform coefficients C(k) and S(k)y respectively, from the encoded signal and generates in 
response to the recovered spectral information a sequence of blocks comprising recovered complex- 
valued frequency-domain transform coefficients ^(fc) of the form ^OcJ-k-yAfk) where 6(k) and A(k) 
30 are obtained from recovered spectral information according to 
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d{k) = cos(^).[e(*) + §{k))] + sin{Jp'[d{J^-k) - S{li-k)l (40) 

A{k) = sin(^).[e(*) + ^(fc))] - cos(:^).ie(^-*) - S(!^.ic)i (41) 

The INDXT and the IMDST constituting the inverse E-TDAC tiansfonn are concurrently 
implemented by an IDFT which generates complex-valued time-domain transform coefficients ^(n) of 
5 the form P(n)+'yf(n) according to 

= jf Ei2(*)«"*^ forO ^ n < V4M (42) 

Tune-domain signal samples x(n) are recovered from the application of the inverse post-transform 
function described above and shown in equations 33 through 36. 

10 

IV. Adaptive-Length Embodiments 
A. Bridge Transform 

As mentioned above, it is desirable tiiat tiie technique which improves transform processing 
efficiency should also permit adaptive selection of Uie transform lengUi. The means and 
15 considaations required to implement an adaptive-transfimn-lengtii co&r are not discussed here, but 
are discussed in International Patent AppUcation PCT/US 91/02512, Publication No. WO 91/16769 
published October 31, 1991). 

Changes in tiie length of either die E-TDAC transform or the 0-TDAC transform may require 
dianges in the phase term m in order to realize time-domain aliasing cancellation. Figure 5 is a 
20 hypotiietical graphical representation of two adjacent overlapped ^-sample length time-domain signal 
sample blocks recovered by an inverse E-TDAC transform, one bloclc recovered from the IMDCT 
and tiie second block recovered from tiie IMDST after syntiiesis windowing but before overlap-add of 
die adjacent blocks has canceUed time-domain aliasing. The representation in tiiis and otiier figures 
does not show individual signal samples, but raflier illustrates only tiie envelope of tiie amplitude of 
25 samples witiiin the windowed signal sample blocks. 

Each recovered signal sample block comprises two components: one component represented in the 
figures by a solid line substantially corresponds to tiie analysis- and syntiiesls-window weighted input 
signal samples, and tiie second component represented in tiie figures by a broken line corresponds to 
tiie analysis- and syntiiesis-window weighted time-domain aliasing distortion. As discussed above, tiie 
30 aliasing component is a time-reversed replica of tiie windowed input signal samples which occurs in 
two separate r^ons. The phase term m for the E-TDAC and tiie 0-TDAC transforms controls tiie 
location of tiie boundary between tiiese two regions. For fixed-lengtii E-TDAC and 0-TDAC 
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transfonns, the boundaiy is located at the mid-point of the signal sample block. The lAase tenn 
required for time^iomain aliasing canceUation under this condition is shown in equation 6. 

Figure 6 is a hypothetical gra^Aical rqiresentation of three time-domain signal sample blocks 
recovaed ftom an inverse E-TDAC transform prior to overlap-add. The first block is an iV-sample 
5 length block wAich has been recovered from the IMDCT. The second and third blocks are 

'AMsample length blocks which have bera recovered from the IMDST. The aliasing component in 
the iV-sample length MDCT block comprises a r^lica of the first half of the signal sample block 
reversed m time about the one-quarter point of the sample block, and a rq)lica of the second half of 
the sampled signal reversed in tune about the three-quarter point of the sample block. If overlap-add 

10 of the second half of the MDCT block and the first MDST ViAT-sample length block shown in 
Figure 6 is to cancel time-domain aliasing, the time-domain aliasing component in the first MDST 
V4iV-sampIe length block must be a replica of the entire ViA^sample length block inverted in sign and 
time-reversed end-for-end. The phase term nt required by the MDST and IMDST transforms to 
produce a time-domain aliasing component with these characteristics is m = l^. 

15 It can be shown that the phase term may be written generally as 

2 

v/hcco ^ is the location of the boundary between time-reversal regions, expressed as the number of 
tune-domain signal samples from the right-hand or trailing edge of the time-domain signal sample 
block. 

20 For example, Figure 7 illustrates two window-weighted time-domain signal sample blocks. The 
right-hand block is UN samples in length. Within this block, the boundary between time-revCTsal 
regions is at a point N/8 samples fcom the right-hand or trailing edge of the block. Thus, the phase 
term m required to cause time-reversal of the aliasing component within each region of the 
UiV-sample block is 

25 m = ±^ — 
2 

W&L this background established, it is now possible to introduce the "bridge transform." A bridge 
transform is a transform Mdiich bridges a shift ftom one transform length to another. For example, as 
shown in Figure 8, suppose the present mvention is called upon to process one block of V^V samples 
followed by another block of UN samples. It is possible to perform a sqiarate transform for each 
30 block. For reasons vdiose explanation exceed the scope of this discussion, a bridge transform 
improves coder performance by instead transforming a single block of ^N samples. 

The bridge transform required to process the %^-sample block shown in Figure 8 may be 
implemented by an FFT to compute the transform for three UN blocks followed by a recombination 
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operatioo* This tedmique is known in the art and is discussed in Oppenheim and Schafer, Digital 
Signal Processing , Englewood Cliffs, NJ.: Prentice-Hall, Inc., 1975, pp. 307-14. The FFT with this 
recombination operation can also be used to concurrently process two E-TDAC bridge transforms in 
the same manner as that briefly discussed above for fixed-length transforms. It is important to note, 
5 however, that concurrent processing in E-TDAC is possible only for a MDCT and MDST which have 
the same length and TDAC phase term. 

For the decoder, the length of the inverse transform may be established by side-information passed 
by the encoder in the encoded signal. The same considerations for ads^tive-length transforms and 
bridge transforms that were discussed above for the forward transform, including the phase term 

10 required for time-domain aliasing cancellation, also apply to the inverse transform. 

The following describes differences between adaptive-length embodiments of the present invention 
and the fixed-length embodiments discussed above. The structure of each adaptive-length 
embodiment is substantially the same as that for a corresponding fixed-length embodiment. The most 
significant differences pertain to the pre- and post-transform functions and to the length and phase 

IS terms of the transform functions. 

In the following discussions, each time-domain signal sample block is defined to be a +2^ samples 
in length, overl^ping tiie immediately prior block by a samples and overlapping the immediately 
subsequent block by b samples. It is assumed that the number of samples in the two overlap intervals 
may vary from block to block. According to the conventions established in the previous discussion, 

20 the bridge transform applied to each time-domain signal sample block is an adaptive-length 
(fl-f-&)-point transform. 

B. E-TDAC Implemented by DCT/DST 
One adaptive-length embodiment corresponds to the fixed-length basic embodiment discussed 
25 above. The forward pre-transform, corresponding to the functions shown in equations 9 and 12 of 
the fixed-length embodiment, generates an alternating sequence of modified sample blocks according 
to 

y(n) = I mod[a+6]) +x([^^-l-n] modla+fc]) | for 0 ^ n < fl^, (44) 

Z{n) = ^ mod[fl+6]) -;c([^l^ -1-n] mod[a+Z^]) | for 0 ^ n < fl^ . (45) 

30 A flowgra^h illustrating this forward pre-transform function for a 16-sample block with a =4 and 
A=12 is shown in Hgure 9. 

The forward transform comprises a DCT and a DST according to 

C(k) = £ y(n)cos(i!^[n+l]) forO ^ * < a+b, (46) 
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^tn*h\ forO ^ ik < a+b. (^7) 



S(.k) = E z(n)siii(^[n+4]) forO ^ * < a+b 



ji-0 



The inverse transform comprises an IDCT and an IDST according to 



Sin) = ^ E ci(k)Cik)cos(^[n^h for 0 ^ n < fl*, (48) 



£(„) - J.g «(]fc)i(*)sm(|^[«-l]) for 0 ^ n < f^. (49) 



5 He inverse post-transform function, corresponding to the functions shown in equations 16, 17, 
19, and 20 for the fixed-length embodunent, recovers time-domain signal samples from time-domain 
transform coefficients according to 

Jt(„) = San-^} mod[«+&]) for 0 ^ n < |, ^ ^ n < a*b, (50) 

m = Sa^-l-n] modla^b]) for f :S n < i^, (51) 

10 Mn) = £([n-^] inod[«+6]) for 0 ^ n < f , ^ n < a^b, (52) 

i(„) = -f ([2^-1-/.] modla*b]) for 1 ^ « < (53) 



C. E-TDAC Implemented by DFT 
Another adaptive-length ^nbodunent corresponds to the fixed-length embodiment of the E-TDAC 
15 transform implemented by a DFT, discussed above. TTie forward pre-transform, corresponding to the 
fonctions shown in equations 21 and 22 for the fixed-length onbodiment, generates an alternating 
sequence of modified sample blocks according to 



(54) 



/,(n) =^|;c([i^+2n]mod[fl+6]) -hx([^-1-2b] modIfl+6]) j for 0^n<^, 

r(n) =^|x([^+2«] mod[a+6]) -x{[?f±-l-2n\ mod[fl+6]) j for 0 ^ « < (55) 

20 The forward E-TDAC transform is implemented by a DFT which generates alternating sets of 
complex-valued ftequency-domain ttansfonn confidents P(k) of the form T(k)+ylJ(k) and R(k) of die 
form V(k)+yW(k) in response to die alternating sequence of modified sample blocks according to 
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m = E PW^"^'*^ for 0 ^ * ^ (56) 

- E 'tH)^?"^'*^ for 0 ^ * ^ fl^. (57) 

The forward post-transform, corresponding to the functions shown in equations 25 and 26 for the 
fixed-length embodiment, generates alternating sets of spectral infonnation according to 

5 CQc) . cos(J* )-7Xfc) + sin(-!i)-J7(*:), (58) 

S(k) = sm(jL) - V(k) - cos(J!L) ' W(k). (59) 

The inverse pre-transfonn function, corresponding to the functions shown in equations 27 through 
30 for the fixed-length embodunent, generates an alternating sequence of blocks comprising recovered 
frequency-domain transform coefficients; one block type comprises recovered complex- valued 
10 coefficients P(k) of the form fl[^)+j •]&(*) and a second block type comprises recovered complex- 
valued coefficients ^(k) of the form 9(k)-^yl^{k), Each frequency-domain transform coefficient is 
obtained according to 

Ilk) = cos(-!* ) . d{k) + siii(J[* ) • C(^-k), (60) 
t){k) = siii( J* ) • t(k) - cosiJlL) • e(f;*-*), (61) 
15 Hk) = cos(-!* ) 'Sif!^-k) * sm(-!* ) -m, (62) 

W(H) - siii(J[*).i(f;*-ife) - cos(-^*).i(*). (63) 

The inverse transform generates an alternating sequence of two types of blocks comprising 
recovered time-domain transform coefficients by applying an IDFT to the alternating sequence of 
frequency-domain transform coefficient blocks; one block type comprises recovered time-domain 
20 transform coefficients p(k) and a second block type comprises recovered time-domain transform 

coefficients r(k). The IDFT used to obtain the recovered time-domain transform coefficients is shown 
in equations 64 and 65; 

Pin) = 2 1> ^(;t)/'*^ for 0 ^ I, < f[J^, (64) 
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a+b ti 2 

The inverse post-transform function, cortesponding to the functions shown in equations 33 through 
36 for the fixed-length embodunent, obtains recovered time-domain signal samples from recovered 
time-domain transform coefficients according to 

5 m = Pil-^^ mod[f^]) = P{[^fi±±^ mod[2±\) for n even, (66) 
2a+&_j_^ 

m = ^([-^-2 3 mod[fI^]) =pa2f±^^ mod[f^]) for n odd. (67) 

_2a+ir 

i(i,) - t{[ ^ ] mod[^]) = /»([2ii:^] mod[^]) for « even, (68) 



i(„) = -/'([-i^^ ] mod[fl*]) = -f([2f!*^] mod[f!I^]) for « odd. 



(69) 



where 0 ^ /» < fl+&. 

10 

D. 0-TDAC Implemented by DST 
The O-TDAC transform utilizes a MDCT of the form 

E{k) = V x(n)cos(2x[Jfc4l-^) forO ^ < (70) 
i:© 2 £1+6 

Efjy = frequency-domain transform coefficient k. 
15 The processing requirements needed to implement this transform can be reduced by a^ilying a 
forward pre-transform function to the time-domain signal samples to generate modified samples e(n)y 
then applying a DST to the modified samples to generate frequency-dommn transform coefficients 
X(k). The forward pre-transform function is 

e{n) =;c([f +n] mod [a+ft]) - x{[l-l-n] mod for 0 ^ n < f , (71) 



20 e(n) = j:([f +n] mod [a+b\) + x([f-l-nl mod \a*b\) for f ^ n < — . (72) 

1 Z Z ^ 



A flowgraph illustrating this forward pre-transform function for a 16-sample block with a=4 and 
6=12 is shown in Figure 10. The minus signs denote terms which are subtractively combined with 
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an associated sample for the fimction shown above in equation 71. 
The forward transform comprises a DST according to 

"2"' 

E{k) = y; tf(/t)sin(iL[*+4][«+4]) forO ^ it < (73) 

The inverse transform comprises an IDST according to 

5 i(n) = -i- E £Wsin(-iL[fc+l]In+i]) forO ^ /i< a+fc (74) 

where i(k) ^ recovered time-domain transform coefficient, and 
E(k) = recovered frequency-domain transform coefRcient k. 

The inverse post-transform function obtains recovered time-domain signal samples x(n) from 
recovered time-domain transform coefficients according to 

10 i(n) = -^([f -1-/I] mod[a+6]) for 0 ^ n < f , (75) 
2 2 

= modlc+ft]) for f ^ n < (76) 

2 2 2 



^(«) = Hli-l-n] inod[a+&]) for ^ » < a+i,. 

2 2 



(T7) 
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CLAIMS 

1. A signal-analysis system for tiie filtering of input samples representing one or more signals 
comprising 

input buffer means for grouping said input samples into time-domain signal sample blocks of 
length ^\^*erein said signal samples are analysis-window weighted samples, and 

analysis means for generating spectral information in response to said time-domain signal 
sample blocks, said spectral information comprising spectral coefficients C(k) and S(k) 
substantially corresponding to the frequency-domain transform coefficients of an Evenly-Stacked 
Time-Domain Aliasing Cancellation transform ^plied to said time-domain signal sample blocks, 
w*erein said spectral coefficients C(k) and S(k) substantially correspond to Modified Discrete 
Cosine Transform coefficients and Modified Discrete Sine Transform coefficients, respectively, 
comprising 

forward pre-transform means for generating modified-sample blocks comprising ^AN 
modified samples by combining one or more pairs of analysis-window weighted samples to 
form said modified samples, and 

forward transform means for generating frequency-domain transform coefficients by 
applying one or more discrete transform functions to said modified-sample blocks. 

2. A signal-analysis system accordmg to claim 1 wherein 

said forward pre-transform means generates first modified-sample blocks comprising modified 

samples y(n) formed firom die additive combination of a pair of analysis-window weighted samples 

x(n) from a respective one of said time-domain signal sample blocks according to 

y(n) «= llx([—-^hl mod N) + xil—-l-n\ mod iV) 1 for 0 ^ n < ^, and said forward 
2 1 4 4 J 2 

pre-transform means generates second modified-sample blocks comprising modified samples z(n) 
formed from the subtractive combination of a pair of analysis-window weighted samples x(n) from 
another respective one of said time-domain signal sample blocks according to 
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z(/f) = 1 mod AO - x([-^-l-n] mod Ao| for 0 ^ /i < and 

said forward transform means generates spectral coefficients C(k) by applying a discrete 
transform function substantially corresponding to a Discrete Cosine Transform function to said 
first modified-sample blocks and generates spectral coefficients S(k) by ^plying a discrete 
transform function substantially corresponding to a Discrete Sine Transform function to said 
second modified-sample blocks. 



3. A signal-analysis system according to claim 1 wherein 

said forward pre-transform means generates first modified-sample blocks comprising modified 
samples p(n) formed from the additive combination of a pair of analysis-window weighted samples 
x(n) from a respective one of said time-domain signal sample blocks according to 
Pin) = l|x([^+2n] mod AO + x([^-l-2n] mod Aoj for 0 ^ n < and said forward 
pre-transform means generates second modified-sample blocks comprising modified samples r(n) 
formed from the subtractive combination of a pair of analysis- window weighted samples x(n) from 
another respective one of said time-domain signal sample blocks according to 
r(n) = i|j:([^+2n] mod AO - x([^-l-2n] mod Ao j for 0 ^ n < 

said forward transfonn means generates a first set of complex-valued frequency-domain 
transform coefficients P(k) of the form Tflc)+yU(k) by applying a discrete transform function 
substantially corresponding to a Discrete Fourier Transform to said first modified-sample blocks 
and generates a second set of complex-valued frequency-domain transform coefficients R(k) of the 
fom V(k)-^yW(k) by applying a discrete transform function substantially corresponding to a 
Discrete Fourier Transform to said second modified-sample blocks, and 
wherein said analysis means further comprises forward post-transform means for generating said 
spectral coefficients C(k) by applymg a forward post-transform function to said first set of complex- 
valued frequency-domain transform coefficients according to C{k) = qos{—) • T\k) + sin(^) • U{k) 
and for generating said spectral coefficients S(k) by applying a forward post-transform function to said 
second set of complex-valued frequency-domain transform coefficients according to 
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S(k) = sm(2p^V(k) - cos(^)-W(*). 
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4. A signal-analysis system according to claim 1 A^erein 

said forward pre-transform means generates said modified-sample blocks comprising complex- 
valued modified samples q(n) of the form p(n)+yr(n) wherein each p(n) is formed from the 
additive combination of a pair of analysis-window weighted samples x(n) from a respective one of 
said time-domain signal sample blocks according to 

p(n) = l|x([^+2n] mod iV) + x([H-l-2n] mod Aoj for 0 ^ n < and each r(n) is 
formed from the subtractive combmation of a pair of analysis-window weighted samples x(n) from 
another respective one of said time-domain signal sample blocks according to 
Tin) « l|x([^+2n] mod iV) - x([M-l-2nl mod ^oj for 0 ^ n < ^, 

said forward transform means generates a set of complex-valued frequency-domain transform 
coefiBcients Q(k) of the form G(k)+yH(k) by applying a discrete transform substantially 
corresponding to a Discrete Fourier Transform to said modified-sample blocks, and 
vvdierein said analysis means further comprises forward post-transform means for generating said 
spectral coefficients C(k) by applying a forward post-transform fonction to sfud set of complex-valued 
frequency-domain transform coefficients according to 

Qk) = iJcos(4j)-[G(Jt) + Gi^'k)] + smA^[H(k) - Hi^-k)] I and for generating said 
2 I iv 2 N 2 J 

spectral coefBcients S(k) by ^plying a forward post-transform fonction to said set of complex-valued 

frequency-domain transform coefficients according to 

Sik) = ^ jcos(^)-[C?(*) - Gi^-m + sm(:^)-[H(Ar) + j. 

5. A signal-analysis system for the filtering of input samples representing one or more signals 
comprising 

input buffer means for grouping said input samples into time-domain signal sample blocks of 
length a+& wherein said length may varies from block to block one or more times, and wherein 
said signal samples are analysis-window weighted samples, 
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forward pre-tiaosfonn means for generating modified-sample blocks comprising ^/i{a-\-b) 
modified samples by combining pairs of analysis-window weighted samples to form said modified 
samples^ and 

forward transform means for generating frequency-domain transform coefficients by applying 
one or more discrete transform functions to said modified-sample blocks. 

6. A signal-analysis system according to claim S wherein 

said forward pre-transform means generates first modified-sample blocks comprising modified 
samples y(n) formed from the additive combination of a pair of analysis-window weighted samples 
x(n) from a respective one of said time-domain signal sample blocks according to 
y{n)^^\^{[?f^^n]mod[a^b])^ for 0 ^ n < and said 

forward pre-transform means generates second modified-sample blocks comprising modified 
samples z(n) formed from the subtractive combination of a pair of analysis-window weighted 
samples x(n) from another respective one of said time-domain signal sample blocks according to 
z(rt)=I jx([^^+n]mod[a+6])-x([Z^-l-n]mod[a+fc^ ^ n < fl^, and 

said forward transform means generates spectral information comprising spectral coefficients 
C(k) and S(k) substantially corresponding to the frequency-domain transform coefficients of an 
Evenly-Stacked Time-Domain Aliasing CanceUation transform applied to said time-domain signal 
sample blocks, >^erein said spectral coefficients C(k) and S(k) substantially correspond to 
Modified IMscrete Cosine Transform coefficients and Modified Discrete Sine Transform 
coefficients, respectively, and wherein said forward transform means generates said spectral 
coefficients C(k) by applying a discrete transform function substantially corresponding to a 
Discrete Cosine Transform function to said first modified-sample blocks and generates said 
spectral coefficients S(k) by applying a discrete transform function substantially corresponding to a 
Discrete Sine Transform function to said second modified-sample blocks. 



7. A signal-analysis system according to claim S wherein 



wo 92/22137 rei/ ua^^/ w/o/ 

.28- 

said forward pre-transfoim means generates first modified-sample blocks comprising modified 
samples p(n) formed from the additive combination of a pair of analysis-window weigjited samples 
x(n) ftom a respective one of said time-domain signal sample blocks according to 
p(«) = ^|A:([^+2n]modtfl+&])+x([2^-l-2n]mod[a+6])|forO ^ n< f^, and said 
forward pre-transform means generates second modified-sample blocks comprising modified 
samples r(n) formed from tfie subtractive combination of a pair of analysis-window weiglited 
samples x(n) ftom anotlier respective one of said time-domain signal sample blocks according to 
r(n) = i |x([^+2«] mod[a+&]) -x([^-l-2n] mod[a+fc]) j for 0 :g « < 

said forward transform means generates a first set of complex-valued frequency-domain 
transform coefficients P(k) of tiie form T(k)-i-yU(k) by applying a discrete transform fimction 
substantially corresponding to a Lfiscrete Fourier Transform to said first modified-sample blocks 
and generates a second set of complex-valued frequency-domain transform coefficients R(k) of the 
form V(k)+yW^) by applying a discrete bsnsform function substantially corresponding to a 
Discr^ Fourier Transform to said second modified-sample blodcs, and 
herein said analysis means furtiier comprises forward post-transform means for generating spectral 
information comprising spectral coefficients C(k) and S(k) substantially corresponding to tiie 
frequency-domain tiansform coefBcients of an Evenly-Stacked Time-Domain AUasing CanceUation 
transform qiplied to said time-domain signal sample blocks, wherein said spectral coefficients C(k) 
and S(k) substantiaUy correspond to Modified Discrete Ctosine Transform coefficients and Modified 
Discrete Sine Transform coefficients, respectively, and wherein said forward transform means 
generates said spectral coefficients C(k) by flying a forward post-transform fanction to said first set 
of complex-valued fiequency-domain transform coefficients according to 



forward post-transform fimction to said second set of complex-valued frequency-domain transform 



C(k) = cos(rlL)'T(k) + sm(-IiL) • U(k) and generates said spectral coefficients S(k) by applying a 



coefficients according to S(k) » sin(-fi) • V(fc) - cos(J^)-W(k). 



8. A signal-analysis system according to claim S wherein 
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said forward pie-transform means gen^ates said modified-sample blocks comprising modified 

samples e(n) formed from the combination of a pair of analysis-window weighted samples x(n) 

from a respective one of said time-domain signal sample blocks according to 

e{n) = +rt] mod - jc([f -l-n] mod [a+ft]) for 0 ^ n < f and 

2 2 2 

e(n) = ;c(lf i-n] mod Ia+6]) + x(lf-l-n] mod for f g n < £1^, and 

2 2 2 2 

said forward transform means generates said spectral information by applying a discrete 
transform function substantially corresponding to a Discrete Sine Transform (DST) function to said 
modified-sample blocks. 



9. A signal-synthesis system for the inverse filtering of spectral information representing one or 
more digital signals comprising 

synthesis means for generating signal samples in response to said spectral information, said 
signal samples substantially corresponding to the time-domain transform coefficients of an Evenly- 
Stacked Time-Domain Aliasing Cancellation transform applied to said spectral information, 
comprising 

inverse pre-transform means for generating sets of frequency-domain transform 
coefficients in response to said spectral information, 

inverse transform means for generating transform blocks comprising time-domain 
transform coefficients by implying an inverse discrete transform function to said sets of 
firequency-domain transform coefficients, and 

inverse post-transform means for generating time-domain signal sample blocks comprising 
N signal samples, wherein one or more pairs of signal samples are generated from a respective 
one of said time-domain transform coefficients, and 

ou^ut means for generating output samples by overlapping pairs of said time-domain signal 
sample blocks and additively combining signal samples from each of said overlapped blocks. 

10. A signal-synthesis system according to claim 9 wherein 
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said inverse transfonn means generates transfonn blocks comprising time^lomain transfonn 
coefficients S(n) by applying an inverse discrete transfonn fanction substantially corresponding to 
an Inverse Discrete Cosine Transfonn function to a respective one of said sets of ftequency- 
domain transfonn coefBcients, and generates transfonn blocks comprising time-domain transform 
5 coefficients z(n) by applying an inverse discrete transfonn fanction substantiaUy corresponding to 

an Inverse Discrete Sine Transfonn fimction to another respective one of said sets of frequency- 
domain transform coefficients, and 

said inverse post-transform means genraates time-domain signal sample blocks comprising 

signal samples x(n) by applying an inverse post-transform fanction to time-domain transfonn 

10 coefficients S(n) according to = Siln-^] mod N) for 0 ^ n< ^ ^ n < N sad 

Ji(n) » ^([1^^-1-n] mod AO for ^ ^ » < ^ applying an inverse post-transfonn 
4 4 4 

fanction to time-domain transform coefficiraits z(n) according to 
m ' £(In-^] mod AO foTO^n<^,^^n< Naod 
ji(n) ' -f ([H-l-n] mod AO for ^ ^ " < 

15 11. A signal-syntiiesis system according to claim 9 wherein said spectral infomation comprises 

blocks of spectral coefficients C:(k) and blocks of spectral coefficients §(k), wherein said spectral 
coefficients t(k) and S(k) substantiaUy correspond to Modified Discrete Cosine Transform coefficients 
and Modified Discrete Sine Transfonn coeffidents of an Evenly-Stacked Tune-Domain Aliasing 
Cancellation transform, respectively, and whodn 
20 said inverse pre-transform means generates first sets of fi«quency-domain transform 

coefficients P{k) of tiiefonn f(k)*ytr{k) according to tm = cos(^)-e(fc) + sm(^) 't{!L-k) 
and t!{k) = sva.i—)-C{k) - cos(—)-di^-k), and generates second sets of frequency-domain 
transfonn coefficients R(k) of die form Vik)+yW(k) according to 
<^(ifc) » cosA'S£-k) * sin(4).5(ife)and ^k) - sinA-i(^-*) - cos(:^).5(*), 
25 said inverse transform means generates transfonn blocks comprising time-domain transform 

coefficients p(n) by flying an inverse discrete transform fimction substantiaUy corresponding to 
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an Inverse Discrete Fourier Transform function to said first sets of frequency-domain transform 
coefficients, and generates transform blocks comprising time-domain transform coefficients f(n) by 
flying an inverse discrete transform function substantially corresponding to an Inverse Discrete 
Fourier Transform function to said second sets of frequency-domain transform coefficients, and 

said inverse post-transform means generates time-domain signal sample blocks comprising 
signal samples x(n) by flying an inverse post-transform function to time-domain transform 
coefficients ^(h; according to Ji{n) = mod[— ]) forO ^ /i < n even and 

6 Z 

^ W = Pil ^^~^ ^^ 1 mod[^]) fotO ^ n < N, n odd, and by applying an inverse post- 

transform function to time-domain transform coefficients r(n) according to 

m = mod[^]) for 0 ^ rt < n even and = -f([l!!L±:^] mod[E]) for 

o 2 8 2 

0 ^ n < iV, rt odd. 



12. A signal-synthesis system according to claim 9 wherein said spectral information comprises 
blocks of spectral coefficients d(k) and blocks of spectral coefficients S(k), wherein said spectral 
coefficients d(k) and S(k) substantially correspond to Modified Discrete Cosine Transform coefficients 
and Modified Discrete Sine Transform coefficients of an Evenly-Stacked Time-Domain Aliasing 
Cancellation transform, respectively, and wherein 

said inverse pre-transfbrm means generates sets of frequency-domain transform coefficients 

Qik) of the form G(k)^yH(k) according to 

Gik) = cosA-[C(*) ^ S(k))] ^ smA'[d(^-k) - 5(^-Jk)]and 
iV yV 2 2 

6(k) = smi^y[d(k) + S(k))] - cos(:^).[d(|-ife) - S(^-kn 

said inverse transform means generates transform blocks comprising complex-valued time- 
domain transform coefficients ^(n) of the form pfnj-hyffn) by applying an inverse discrete 
transform substantially corresponding to an Inverse Discrete Fourier Transform, and 

said inverse post-transform means generates time-domain signal sample blocks comprising 

signal samples x(n) by applying an inverse post-transform function to p(n) according to 

m = pil^!^] mod[:J]) for 0 < n < M n even and Jt{n) = pa^I^ltl^] mod[E]) for 
o 2 8 2 
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0 ^ n < //, n odd, and by flying an inverse post-transfonn function to r(n) according to 
Jtin) = m ^"'^^ ] mod[~]) for 0 ^ n < M /« even and m - -/>([ 3N-4-4n j ^^j^j^ 
0 ^ n < N,n odd, 

13. A signal-synthesis system for the inverse filtering of spectral information rq)resenting one or 
more digital signals comprising 

inverse pre-transform means for generating sets of ftequency-domain transform coeflBcients in 
response to said spectral information, said sets comprising coefficients wherein the 

number of coefficients may varies from set to set one or more times, 

inverse transform means for generating transform blocks comprising time-domain transform 
coefRcimts by applying one or more inverse discrete transform functions to said sets of frequency- 
domain transform coefBcirats, 

inverse post-transform means for generating time-domain signal sample blocks comprising 
a+i signal samples, M*erein a pair of signal samples is generated from a respective one of said 
time-domain transform coefficients, and 

ou^ut means for generating ou^ut samples by overlaying pairs of said time-domain signal 
sample blocks and additively combining signal samples fipom each of said overlapped blocle, 

14. A signal-synthesis system according to claim 13 ^^erein 

said inverse transform means generates transform blocks comprising time-domain transform 
coefficients by applying an inverse discrete transform function substantially corresponding to 
an Inverse Discrete Cosine Transform function to a respective one of said sets of frequency- 
domain transform coefficients, and transform blocks comprising time-domain transform 
coefficients z(n) by applying an inverse discrete transform function substantially corresponding to 
an Inverse Discrete Sine Transform function to another respective one of said sets of frequency- 
domain transform coefficients, and 

said inverse post-transform means generates time-domain signal sample blocks comprising 
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signal samples i(n) by applying an mverse post-tiansfonn function to time-domain ttansfonn 
coefficients according to 

m = S(ln-^^1 mod[a+6]) for 0 ^ n < ^, ^ n < a*b and 

^(«) = ^d-^^^-l-n] mod[a+ft]) for f ^ n < and by applying an inverse post- 

2 2 2 

transform function to time-domain transform coefficients z(n) according to 

f («) - f([n--?i^] mod[a+i]) for 0 ^ n < f , 1^ ^ n < a+b and 
2 2 2 

m = mod[a+6]) for f ^ n < 

2 2 2 



IS. A signal-synthesis system according to claim 13 wherein said spectral information comprises 
blocks of spectral coefficients d(k) and blocks of spectral coefficients S(k), wherein said spectral 
coefficients C^) and S(k) substantially correspond to Modified Discrete Cosine Transform coefficients 
and Modified I^rete Sine Transform coefficients of an Evenly-Stacked Time-Domain Aliasing 
Cancellation transform, respectively, and wherein 

said inverse pre-transform means generates first sets of frequency-domain transform 

coefficients P(k) of the form f(Jfc)+j'I/(ifc) according to 

TXk) = cos(-^* )-C(Jfc) + sin(-!* and 
a-^b a-^b 2 

t^(k) = sin(-I— ) • C^k) - cos(-l!L) • and generates second sets of frequency-domain 

a+6 a +6 2 

transform coefficients j^(k) of the form ^ik)'^yW{k) according to 

V(k) = cosiJL)^S(—-k) + sin(J[i)-^(it) and 
a-^b 2 a+b 

W(k) = sin(J*).5(f;*-*) - cos(_2* 
a+b 2 a+b 

said inverse transform means generates transform blocks comprising time-domain transform 
coefficients p(n) by applying an inverse discrete transform function substantially corresponding to 
an Inverse Discrete Fourier Transform function to said first sets of frequency-domain transform 
coefficients, and generates transform blocks comprising time-domain transform coefficients f(n) by 
^plying an inverse discrete transform function substantially corresponding to an Inverse Discrete 
Fourier Transform function to said second sets of frequency-domain transform coefficients, and 

said inverse post-transform means generates time-domain signal sample blocks comprising 
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sdgnal samples x(n) by applying an inverse post-transfonn fiincUon to time-domain transfonn 

coefficients p(n) accoiding to i («) = pjl^'^'^] mod[^]) for 0 ^ « < a+b, n even and 

i(„) « ^([if±±^] modlfl*]) for 0 S » < o+i, n odd, and by applying an inverse post- 

ttansfonn function to time-domain transfonn coefficients f(n) according to 

jt{n) » /*([^"^"^] mod[fJ^]) fbrO ^ n < a+Z^, n even and 
4 2 

je(„) = -/.([2g-^&-2-2Hj „od[f!;^]) for 0 ^ n< a+b, n odd. 
4 2 

16. A signal-synthesis system according to claim 13 wherein 

said inverse transform means generates tune-domain transform coefficients S(n) by ^plying an 
inverse discrete transform function substantially corresponding to an Inverse Discrete Sine Transform 
to said sets of frequency-domain transfonn coefficients, and 

\^*erein said inverse post-transform means generates time-domain signal sample blocks comprising 

signal samples x(n) by flying an inverse post-transform function to time-domain transfonn 

coefficients i(n) according to je(n) « -«([|-l-n] mod[a+6]) for 0 ^ n < f , 

=2([n--J] mod[a+6]) for ^ ^ n < and 
2 2 2 

=g([£-l-/i]mod[a+6]) for 3^ ^ n < a+6. 

2 2 
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